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Abstract. We present estimates of the amplitudes of in- 
trinsically stable stochastically excited radial oscillations 
in stars near the main sequence. The amplitudes are de- 
termined by the balance between acoustical energy gener- 
ation by turbulent convection (the Lighthill mechanism) 
and linear damping. Convection is treated with a time- 
dependent, nonlocal, mixing-length model, which includes 
both convective heat flux and turbulent pressure in both 
the equilibrium model and the pulsations. Velocity and lu- 
minosity amplitudes are computed for stars with masses 
between 0.9 Mq and 2.0 in the vicinity of the main 
sequence, for various metallicities and convection param- 
eters. As in previous studies, the amplitudes are found 
to increase with stellar mass, and therefore with luminos- 
ity. Amongst those stars that are pulsationally stable, the 
largest amplitudes are predicted for a 1.6 model of 
spectral type F2; the values are approximately 15 times 
larger than those measured in the Sun. 
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1. Introduction 

The stability of solar-like p modes depends mainly on the 
interaction of the oscillations with radiation and convec- 
tion in the outer envelope. The most plausible explana- 
tions for the occurrence of such oscillations are either in- 
trinsic thermal overstability or stochastic excitation of sta- 
ble modes by turbulent convection. Whatever the mecha- 
nism, the energy flow from radiation and convection into 
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and out of the p modes takes place very near the surface 
(e.g. Goode, Cough & Kosovichev 1992). Sun-hke stars 
possess surface convection zones, and it is in these zones, 
where the energy is transported principally by the turbu- 
lence, that most of the driving takes place. Mode stability 
is governed not only by the perturbations in the radiative 
fluxes (via the K-mechanism) but also by the perturba- 
tions in the turbulent fluxes (heat and momentum). The 
study of mode stability therefore demands a theory for 
convection that includes the interaction of the turbulent 
velocity field with the pulsation. 

Thermal overstability of pulsations arises from an ex- 
change of energy between the oscillatory motion of the 
stellar matter, the turbulence and the radiation field. Such 
overstability has been suggested as a possible mechanism 
for the excitation of solar oscillations by Ulrich (1970a) 
and Antia, Chitre & Cough (1988). If solar p modes were 
indeed overstable, some nonlinear mechanism must limit 
their amplitudes to the values that are observed. Non- 
linear coupling to other, stable modes was considered by 
Kumar & Goldreich (1989), who estimated that the en- 
ergy drain through three-mode coupling would occur at a 
rate too low to extract the energy gained from the over- 
stability at the appropriate amplitude. Similar estimates 
of nonlinear self limiting are also too weak. The satura- 
tion of mode amplitudes at the observed levels therefore 
remains a mystery if overstability provides the origin of 
solar pulsations. 

The problem of identifying a saturation mechanism 
does not arise if the modes are intrinsically stable. Such 
modes can be stochastically excited by the turbulent con- 
vection. The process can be regarded as multipole acous- 
tical radiation (e.g. Unno 1964). For solar- like stars, the 
acoustic noise generated by convection in the star's reso- 
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nant cavity may be manifest as an ensemble of p modes 
over a wide band in frequency (Goldreich & Keeley 1977b). 
The amplitudes are determined by the balance between 
the excitation and damping, and are expected to be rather 
low. The turbulent-excitation model predicts not only 
the right order of magnitude for the p-mode amplitudes 
(Gough 1980), but it also explains the observation that 
millions of modes arc excited simultaneously. Moreover, 
recent observations (Toutain & Frohlich 1992, Goode & 
Strous 1996; Chaplin et al. 1998) also support a stochas- 
tic origin. This second explanation seems therefore to be 
the more likely, and is the one that we shall adopt here. 

To date, Christensen-Dalsgaard & Prandsen (1983b) 
have made the only predictions of amplitudes of solar-like 
oscillations in other stars. They obtained amplitudes of 
modes by postulating equipartition between the energy of 
an oscillation mode and the kinetic energy in one convec- 
tive eddy having the same turnover time as the period 
of the oscillation. This simple formula for excitation was 
proposed by Goldreich & Keeley (1977b), who used it to 
estimate amplitudes for the solar case, assuming damp- 
ing rates determined solely by a scalar turbulent viscos- 
ity (Goldreich & Keeley 1977a). In the calculations of 
Christensen-Dalsgaard & Frandsen (1983b), however, ra- 
dial eigenfunctions of model envelopes were computed by 
solving the equations of linear nonadiabatic oscillation, 
although they also neglected turbulent pressure and set 
the Lagrangian perturbation to the convective heat flux 
to zero. They found velocity and luminosity amplitudes 
to increase with age, and with increasing mass along the 
main sequence. 

Balmforth (1992a) improved the calculation by intro- 
ducing Cough's (1976, 1977) nonlocal, time-dependent 
mixing-length model for convection, using the Eddington 
approximation to radiative transfer for both the equilib- 
rium structure and the pulsations. He calculated damping 
rates for the solar case and found all modes to be stable. 
Here we continue Balmforth's investigation and study the 
oscillations of main-sequence stars, delimiting the region 
in the HR diagram for stars with stable modes. Prelim- 
inary results of the calculations have been presented by 
Houdek et al. (1995). 

According to Libbrecht et al. (1986) the observed oscil- 
lation properties of low-degree modes depend little on the 
value of I. This is to be expected because the excitation 
and damping mechanisms are significant only very close 
to the surface, where the vertical scale is much less than 
the horizontal scale of oscillations and when I is low the 
modal inertia is quite insensitive to degree L It is therefore 
adequate to simplify the calculations, by analysing only 
radial modes of oscillation. The results are applicable to 
all modes of moderately low degree. 



2. Observational projects 

Observations of oscillation properties of stars other than 
the Sun provide important information for testing the the- 
ory of stellar evolution. A critical problem with the detec- 
tion of oscillations in solar-type stars, however, is their 
very small amplitudes, of the order of 1 ms~^ or less. For 
the Sun the observed velocity variations in disc-integrated 
light have values < 20cms~^ (e.g. Grec et al. 1983; Lib- 
brecht & Woodard 1991; Chaplin et al. 1998). To detect 
similar variations in distant stars is therefore a challenging 
task, requiring observations to be made with the utmost 
precision. 

Three observing techniques for detecting such oscilla- 
tions have been developed so far. The first is to search for 
periodic Doppler shifts of spectral lines (e.g. Kennelley 

1995) . However, the most successful result by this method 
has been the determination of only an upper bound to 
oscillation amplitudes in some of the brightest stars (e.g. 
Brown & Gilliland 1990; Brown et al. 1991; Mosser et al. 
1998). The second method is to look for periodic bright- 
ness fluctuations using photometry. When used with area 
detectors such as CCDs, this method has a clear advan- 
tage over spectroscopic techniques because it permits one 
to observe large ensembles of stars simultaneously (e.g. 
Gilliland 1995). Using differential CCD photometry with 
seven 4m-class telescopes, Gilliland et al. (1993) obtained 
upper limits of luminosity amplitudes of possible oscilla- 
tions in twelve stars in M67. The third method, introduced 
by Kjeldsen & Bedding (1995), measures temperature fluc- 
tuations induced by stellar oscillations via their effects on 
the equivalent width of the Balmer lines (Bedding et al. 

1996) . This technique has yielded a possible detection of 
solar-like oscillations in the sub-giant rj Boo (Kjeldsen et 
al. 1995). Although it is currently restricted to observing 
isolated stars, the equivalent-width method is insensitive 
to atmospheric scintillation, and attains a substantially 
better signal-to-noise ratio than do the other two ground- 
based methods. 

The limitations of ground-based observational tech- 
niques in asteroseismology have been addressed by Frand- 
sen (1992) and Gilhland (1995), both of whom argued that 
seismology can be applied to distant solar-type stars only 
by observing them from space. The elimination of atmo- 
spheric noise and the possibility of obtaining long con- 
tinuous data sequences will provide information of much 
higher quality than from any ground-based method. Sev- 
eral asteroseismological space projects are in preparation, 
such as the French project COROT^ (Catala et al. 1995), 
the Danish project MONS^ (Kjeldsen & Bedding 1998) 
and the Canadian project MOST^ (Matthews 1998). 
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When preparing an observing campaign, it is helpful 
in the selection of target stars to have a good prediction 
of the amplitude of the signal one will be trying to ob- 
serve. Measurements of mode lifetimes (damping rates) 
and the variation of oscillation amplitudes, which dc^pcnd, 
in part, on stellar parameters, provide invaluable insights 
into the mechanisms that excite solar-like oscillations. It 
is hoped that future observations will provide these cru- 
cial measurements. The aim of this paper is to provide a 
systematic survey of the oscillation properties in view of 
these upcoming observational projects. 

3. Time-dependent convection 

In order to describe the turbulent fluxes in a time- 
dependent envelope, various phenomenological mixing- 
length models have been proposed (for a detailed dis- 
cussion see Balmforth 1992a; Houdek 1996). In envelopes 
that do not pulsate, the various guises of local mixing- 
length method are essentially similar, once their intrin- 
sic parameters have been calibrated (their formulations 
may be interpreted as providing interpolation formulae 
between the two limits of efficient and inefficient convec- 
tion; Gough & Weiss 1976). This is not so of formulations 
of time-dependent mixing-length prescriptions, in which 
the details of the phenomenological model inffuence the 
predictions of stability (Balmforth 1992a). Moreover, lo- 
cal theories are plagued by some fundamental inconsisten- 
cies, which we shall describe presently. For these reasons, 
we prefer to use a nonlocal version (Gough 1976) which is 
summarized below. 

3.1. Local mixing-length model 

Local mixing-length models (e.g. Bohm-Vitense 1958; Ul- 
rich 1970b; Gough 1977) still provide the almost univer- 
sal method for computing the stratification of convection 
zones in stellar models. One of the major drawbacks of a 
local approach is the assumption that the characteristic 
length scale I must be shorter than any scale associated 
with the structure of the star. This condition is certainly 
violated in solar-like stars and in red giants, where cali- 
brated evolution calculations yield a typical value for the 
mixing- length parameter a = IjHp that is of order unity, 
where is the pressure scale height. (Since a is normally 
thought to be essentially invariable, it follows that the con- 
dition must be thought to be violated in all stars.) This 
implies that fluid properties vary significantly over the ex- 
tent of a convective element; the super adiabatic gradient 
can vary on a scale much shorter than i. 

In many computations of stellar envelope models, 
the turbulent pressure pt associated with the Reynolds 
stresses has been ignored. However, several investigations 
(Baker & Gough 1979; Rosenthal et al. 1995; Canuto & 
Christensen-Dalsgaard 1998) suggest that the momentum 



flux provides a substantial fraction of the hydrostatic sup- 
port in the equilibrium model, and should therefore not 
be neglected. In a local model, inclusion of the momen- 
tum flux pt = pw^ {p being the mean density and w the 
rms vertical velocity of the convective elements) leads to 
singular points in the equation for the turbulent pressure 
gradient at the edges of the convective regions (e.g. Gough 

1976) . This issue demands careful consideration, although 
it can be circumvented by adopting an additional approxi- 
mation to reduce the order of the equations of hydrostatic 
support, which are then not strictly consistent with the 
formulation of the theory (Henyey, Vardya & Bodenheimer 
1965; Baker & Gough 1979). 

In a time-dependent description, the details of the phe- 
nomenological model become important to the linear sta- 
bility calculation; the nonlocal model that we employ in 
our calculations is one that is based on Gough's (1965, 

1977) local mixing-length model. In that formulation, the 
turbulent eddies that support the heat and momentum 
fluxes evolve in a pulsating environment. Explicit con- 
sideration is given to the phase of pulsation at the birth 
of each transitorily coherent eddy, and to how the eddy 
adjusts to the temporally varying environment. Other 
mixing-length models (e.g. Unno 1967, Xiong 1989) em- 
phasize other aspects of the dynamics, which presumably 
influence differently the pulsational stability (Balmforth 
1992a). 

Another major drawback of local theory is that it fails 
to treat the convective dynamics across extended eddies 
in a physically plausible fashion: in deeper parts of the 
convection zone, where the stratiflcation is almost adia- 
batic, convective heat transport is very efficacious: radia- 
tive diffusion is unimportant, and the perturbation of the 
heat flux is dominated by advection of temperature fluc- 
tuations. In this limit, the temperature fluctuation can be 
described approximately by a diffusion equation in which 
the diffusivity is imaginary (Baker & Gough 1979; Gonczi 
& Osaki 1980). Rapid spatial oscillations of the eigenfimc- 
tions result. The problem this introduces can at best be 
thought of as one of numerical resolution, which is par- 
ticularly severe in layers where the stratification is very 
close to being adiabatic. But worse, it signifies a complete 
failure of the model, since the wavelength of the spatial 
oscillation is very much smaller than the mixing length 
(which is supposedly the smallest lengthscale permitted 
by the model). 

3.2. Nonlocal mixing-length model 

The obvious drawbacks of a local formulation of the 
mixing-length approach can be removed by an appropri- 
ate nonlocal generalization. Account can be taken of the 
finite size of a convective eddy by averaging spatially the 
representative value of a physical variable throughout the 
eddy. Spiegel (1963), for example, proposed a nonlocal de- 
scription based on the concept of an eddy phase space, and 
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derived an equation for the convective flux which is similar 
to a radiative-transfer equation. The solution of this trans- 
fer equation yields an integral expression that converts 
the usual ordinary diflferential equations describing stel- 
lar models into integro-differential equations. The solution 
to the transfer equation can be approximated by taking 
moments, closing the hierarchy at second order with the 
Eddington approximation (Gough 1976). In this approx- 
imation, the nonlocal generalization of the mixing-length 
formulation is governed by three second-order differential 
equations for the spatially averaged convective fluxes Fc 
and pt and for the average superadiabatic temperature 
gradient experienced by an eddy. 

The procedure introduces two more parameters, a and 
b, which characterize respectively the spatial coherence of 
the ensemble of eddies contributing to the total heat and 
momentum fluxes, and the extent over which the turbu- 
lent eddies experience an average of the local stratification. 
Theory suggests approximate values for these parameters, 
but it is arguably better to treat them as free. Roughly 
speaking, the parameters control the degree of "nonlocal- 
ity" of convection; low values imply highly nonlocal solu- 
tions, and in the limit a, 5 ^ oo the system of equations 
reduces to the local formulation (except near the bound- 
aries of the convection zone, where the equations are sin- 
gular). Balmforth (1992a) investigated the effect of the 
parameters a and b on the turbulent fluxes in the solar 
case, and Tooth & Gough (1989) tried to calibrate a and 
b by comparing theory with laboratory experiments. 

4. Model computations 

The basic model calculations reported in this paper are as 
described by Balmforth (1992a). In particular, we incorpo- 
rate turbulent pressure in the equilibrium model envelope. 
The computation proc;eeds by iteration, from a trial solu- 
tion obtained by integrating inwards from an optical depth 
of T = 10~* and ending at a radius fraction 0.2, using lo- 
c;al mixing-length theory and the diffusion approximation 
to radiative transfer; the approximation to the turbulent 
pressure used by Baker & Gough (1979) is adopted to ob- 
viate singular points at the edges of the convection zones. 
The entire envelope is then re-integrated using the equa- 
tions appropriate to the nonlocal mixing-length theory, 
using the Eddington approximation to radiative trans- 
fer (Unno & Spiegel 1966). The atmosphere is treated 
as being grey, and is assumed to be plane parallel. The 
temperature gradient is corrected by using a r-dependent 
varying Eddington factor (Auer & Mihalas 1970) derived 
from model C of Vernazza, Avrett & Loeser (1981) in the 
manner of Balmforth (1992a). Opacities were obtained 
from the latest OPAL tables (Iglesias & Rogers 1996), 
supplemented at low temperature by tables from Kurucz 
(1991). Interpolation in these tables was carried out using 
birational splines (Houdek & Rogl 1996). The equation 



of state included a detailed treatment of the ionization 
of C, N, and O, and a treatment of the first ionization 
of the next seven most abundant elements (Christensen- 
Dalsgaard 1982), as well as 'pressure ionization' by the 
method of Egglcton, Faulkner & Flannery (1973); elec- 
trons were treated with relativistic Fermi-Dirac statistics. 
In the pulsation model the boundary conditions used are 
essentially those of Baker & Kippenhahn (1965), but sup- 
plemented by appropriate conditions on the variables of 
the nonlocal mixing- length theory (Balmforth 1992a). The 
outer boundary conditions were applied at the tempera- 
ture minimum, the mechanical condition being consistent 
with a perfectly reflecting surface; at the base of the en- 
velope, conditions of adiabaticity and vanishing displace- 
ment were imposed. 

The linearized pulsation equations were solved with a 
second-order accuracy Newton-Raphson-Kantorovich al- 
gorithm (Baker, Moore & Spiegel 1971; Gough, Spiegel 
& Toomre 1974). With this algorithm the eigenfunctions 
and eigenvalues can be computed simultaneously; how- 
ever, one has to provide a proper trial solution. This can 
be obtained by solving flrst the adiabatic pulsation equa- 
tions and then applying a quasi-adiabatic approximation* 
to complete the nonadiabatic system. The detailed equa- 
tions describing the equilibrium and pulsation model have 
been discussed by Balmforth (1992a) and Houdek (1996). 

5. Damping rates 

Were solar p modes to be genuinely linear and stable, 
their power spectrum could be described in terms of an 
ensemble of intrinsically damped, stochastically driven, 
simple-harmonic oscillators (Batchelor 1956; Christensen- 
Dalsgaard, Gough & Libbrecht 1989) provided the back- 
ground equilibrium state of the star is independent of 
time; if we further assume that mode phase fluctuations 
contribute negligibly to the width of the spectral lines, the 
intrinsic damping rates of the modes could then be deter- 
mined observationally from measurements of the pulsation 
linewidths. The linewidths are obtained, to a first approxi- 
mation, by fitting Lorentzian profile functions to the spec- 
tral peaks of the observed power spectrum. Higher approx- 
imations demand a more detailed description of the excita- 
tion and damping (cf. Jefferies et al. 1988; Gabriel 1998). 
Continuous observations over many periods are required. 
Observations from GONG (Harvey et al. 1996), and from 
the SOI-MDI (Scherrer, Hoeksema & Bush 1991), VIRGO 
(Frohhch et al. 1995) and GOLF (Gabriel et al. 1991) in- 
struments aboard the SOHO spacecraft have already pro- 
vided such high-quality data for the Sun. 

It should be noted that the observed power spectrum 
of solar pulsations is complicated by the beating of closely 

Quasi-adiabatic approximations adopt adiabatic eigen- 
functions for evaluating the thermal variables for comput- 
ing mode stability. 
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incoherent 
scattering; prin- 
cipally at the 
horizontally in- 
homogeneous 
superadiabatic 
boundary layer 
(Goldreich & 
Murray, 1994) 



transmission of 
waves into the 
atmosphere 

(Balmforth & 
Gough, 1990b) 



modulation 
of the turbu- 
lent momen- 
tum flux P( 
by the pulsa- 
tion (Gough, 
1977) 



radiative damping 
due to nonadiabatic 
effects, e.g. K-mech., 
departure from radia- 
tive equilibrium in 
the mean state 
(Christensen- 
Dalsgaard & Frandsen, 
1983a) 



modulation of 
the turbulent 
heat flux by 
the pulsation 
(Gough, 1977) 



Fig. 1. Physical processes contributing to the 
linear damping rate r]. They can be associated 

with the effects arising from the momentum bal- 
ance (77dyn) and from the thermal energy balance 
(rig). The contributions r^scatt and r/icak arc in 
parentheses because they have not been taken 
into account in the computations reported in this 
paper. The influence of Reynolds stresses on so- 
lar modes, contributing to r]t , has been treated by 
Goldreich & Keeley (1977a) in the manner of a 
time-independent scalar turbulent viscosity. The 
width of the line in the Fourier power spectrum 
of the oscillations is influenced also by nonlin- 
earities, both those coupling a mode to others 
(Kumar & Goldreich, 1989) and those intrinsic 
to the mode itself. 



spaced modes (e.g. Christensen-Dalsgaard & Gough 1982; 
Hill et al. 1996) and by amplitude variations and phase 
wandering resulting from the interaction with the tur- 
bulent convection. These effects might also be responsi- 
ble for the observed asymmetries in the p-mode line pro- 
files. The study of how to disentangle the manifestations 
of these phenomena from manifestations of the intrinsic 
mode parameters is only in its infancy (Chang, Gough 
k Sekii 1997; Roxburgh & Vorontsov 1997; Gabriel 1998; 
Nigain ct al. 1998; Rast & Bogdan 1998). It is important 
for analysing not only the oscillations of the Sun, but the 
oscillations of any star with a rich spectrum of frequencies. 

5.1. Processes contributing to intrinsic linewidths 

Basically, the damping of stellar oscillations arises from 
two sources: processes influencing the momentum balance, 
and processes influencing the thermal energy equation. 
Each of these contributions can be divided further ac- 
cording to their physical origin, as illustrated in Fig. 1. 
A detailed discussion of the processes has been given by 
Houdek (1996). Here wc limit the discussion to those that 
are modelled in our computations. 

Nonadiabatic radiative processes can contribute to 
both the driving and the damping of the pulsations. In 
solar-type stars the zones of ionization lie well inside the 
regions of efScacious convection, and the conventional k- 
mechanism provides only a relatively small contribution to 
the driving. Radiative damping in the atmosphere is not 
necessarily small, and requires a more accurate treatment 
of radiative transfer than the diffusion approximation. 
Christensen-Dalsgaard & Frandsen (1983a) have shown 
that the use of the grey Eddington approximation, when 
applied correctly, does not introduce too large an error in 
the calculation of the damping rates. Ftirthermore, they 
have demonstrated that for stability calculations depar- 
tures from radiative equilibrium in the mean state must 
not be neglected: in the upper boundary layer of the con- 
vection zone, where there is a transition from convective 
to radiative energy transport, radiative equilibrium is no 
longer maintained. Thus the mean intensity J is not equal 



to the Planck function B. In particular, by perturbing 
the equations describing the radiation field in the Edding- 
ton approximation, one obtains (Christensen-Dalsgaard & 
Frandsen 1983a) 

(5(-divFr) = ^ttk[6B - 6J+—{Bo - Jq)] , (1) 

p K 

where p and k denote the density and opacity, respectively, 
the operator 6 denotes a Lagrangian perturbation, and 
the subscript denotes an equilibrium quantity. The last 
term in Eq. (1) describes the departure from radiative 
equilibrium in the mean state; it is not everywhere small, 
yet it has been ignored in most stability calculations so 
far. 

Vibrational stability is influenced further by the ex- 
change of energy between the pulsation and the turbulent 
velocity field. The exchange arises cither via the pulsation- 
ally perturbed convective heat fiux, or directly through 
dynamical effects of the fluctuating Reynolds stresses. In 
fact, it is the modulation of the turbulent fluxes by the pul- 
sations that seems to be the predominant mechanism re- 
sponsible for the driving and damping of solar-type acous- 
tic modes. 

Nonadiabatic processes attributed to the modulation 
of the convective heat flux by the pulsation are accounted 

for by the contribution 77conv to the total damping rates 
(see Fig. 1). This contribution is related to the way that 
convection modulates large-scale temperature perturba- 
tions induced by the pulsations, which influences pulsa- 
tional stability substantially. The manner in which it does 
so, together with the conventional K-mechanism, is dis- 
cussed by Balmforth (1992a). It appears to have a signif- 
icant destabilizing influence on the pulsations (Balmforth 
& Gough 1990a). 

It was flrst reported by Gough (1980) that the dynam- 
ical effects arising from the turbulent momentum flux per- 
ttirbations 5pi contribute significantly to the damping r]t ■ 
Detailed analyses (Balmforth 1992a) reveal how damping 
is controlled largely by the phase difference between the 
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Fig. 2. Linear damping rates r) for the Sun as function of fre- 
quency. The values chosen for the convection parameters are 
a = 1.8 and = \? = 600. 

turbulent pressure perturbation 3i\ and the density per- 
turbation bp. Turbulent pressure fluctuations must not, 
therefore, be neglected in stability analyses of solar-type 
p modes. 

The results presentcid here were obtained from com- 
putations including the physics describing r^rad, ''7conv and 
?7t. The nonadiabatic contributions ?7rad and ryconv may be 
associated with the thermodynamics of the gas, and ac- 
cordingly we couple them into = rjrad + Vconv (see also 
Fig. 1). 

5.2. Theoretical damping rates 

Damping rates are computed as the imaginary part uii of 
the complex eigenfrequency ui = Ur + itOj, obtained from 
solving the fully nonadiabatic pulsation equations. Balm- 
forth (1992a) computed damping rates for the Sun, and 
reported that he found all modes to be stable, with damp- 
ing rates agreeing tolerably with observation for frequen- 
cies between 2 mHz and 4 mHz. Below and above this 
frequency range the theoretical damping rates are smaller 
than observations would suggest. Damping arising from 
incoherent scattering r^scatt (Goldreich & Murray 1994, see 
Fig. 1), which may remove the discrepancy both at low and 
at high frequencies, is not modelled in our calculations. 

Figure 2 displays the damping rates and their contri- 
butions arising from the gas and turbulent pressure per- 
turbations for a solar envelope model. Damping is much 
augmented by the turbulent pressure perturbation 6pt; it 
is only at the highest frequencies that the nonadiabatic 
contribution rjg to damping of solar p modes exceeds that 
from the turbulent pressure 

The total damping rate r] (solid curve), plotted as a 
function of cyclic frequency i> = lu,;/2'k, is characteristi- 
cally flat at frequencies near 2.8 mHz (see Fig. 2). This fea- 
ture is also observed in solar linewidth measurements (e.g. 
Libbrecht 1988, Appourchaux et al. 1998, Chaplin et al. 



I I 

9.72 Gy 




1 2 3 4 5 6 

1/ [mHz] 



Fig. 3. Damping rates for an evolving 1 Mq star as function of 
frequency. The results are displayed for models with ages=(0, 
2.49, 3.96, 4.55, 6.19, 7.00, 8.03, 9.02, 9.72) Gy. The thick curve 
indicates the results for the present Sun. Values a = 1.8, = 
= 600 for the convection parameters have been used. 
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Fig. 4. Damping rates for an evolving 1.45 Mq star as a func- 
tion of frequency. The results are depicted for models with 
ages=(0, 0.96, 1.38, 1.72, 2.00, 2.44) Gy. Values a = 2.0, = 
900, = 2000 for the convection parameters have been used. 

1998). At these frequencies the net damping is reduced 
particularly by radiative processes in the upper superadi- 
abatic boundary layer of the convection zone, which are 
locally destabilizing. 

The damping rates for an evolving 1 Mq star are de- 
picted in Fig. 3. Damping rates generally increase with in- 
creasing age, particularly for low- and high-order modes. 
For modes of intermediate order the flattening of the 
damping-rate curve becomes more pronounced as the star 
evolves, and turns into a locally concave function at about 
the solar age. The maximum value of the superadiabatic 
temperature gradient of a 1 Mq star increases by approx- 
imately 24% along the main sequence, promoting the de- 
pression in the damping rates. 
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Fig. 5. Damping rates for ZAMS models as functions of fre- 
quency. The results are displayed for models with M=(0.95, 
1.00, 1.05, 1.10, 1.15, 1.20, 1.25, 1.30, 1.35, 1.40, 1.45, 2.00) 
Mq. For the convection parameters the values a = 2.0, = 
900, Ip = 2000 have been used. 

A similar behaviour of the damping rates is obtained 
for more massive stars, as indicated for the evolving 
1.45 M0 star depicted in Fig. 4. For these more massive 
stars, larger values of the nonlocality parameters, a? = 900 
and 6^ = 2000, were adopted. This selection ensured that 
all the radial modes were stable, in line with our work- 
ing hypothesis that stochastic excitation underlies the ap- 
pearance of solar-like oscillations (but see section 9). The 
depression in the damping rates is more pronounced for 
these stars, even at the ZAMS. This trend may be seen 
even more obviously in Fig. 5, where damping rates are 
depicted for stars with increasing mass along the ZAMS. 

The functional dependence of rj on stellar parameters 
was determined approximately by Goldreich & Kumar 
(1991); they derived an ordcr-of- magnitude estimate for 
the damping rates accounting roughly for the effects of 
radiative damping and convective dynamics. The radial 
modes were treated in the polytropic approximation to 
the outer layers and convection was described by standard 
(unperturbed) local mixing-length theory (Bohm-Vitense 
1958). They obtained the expression 



(2) 



where L is the luminosity and Cg denotes the adiabatic 
sound speed at the photosphere (which we define at the 
level where the temperature is equal to the effective tem- 
perature), lo, is the mode inertia, and uj^ is the acousti- 
cal cut-off frequency in an isothermal atmosphere (Lamb 
1909), 



Wc = 



2H„ 



(3) 



where c denotes the sound speed. The inertia is usually 
defined such that it represents the coefficient of propor- 



tionality between the energy in the mode and the square of 
the velocity amplitude of the associated disturbance in the 
surface layers of the star. The surface of the star, however, 
commonly lies in a region where the mode is evanescent, 
and in that case 1^ is more usefully regarded as a measure 
of evanescence, representing a property of the eigenfunc- 
tion above the upper turning point (Gough 1995). For the 
case of a linear adiabatic mode of stellar oscillation, which 
can be represented by an undamped harmonic oscillator, 
the mode inertia 1^ can be defined in terms of the total 
(kinetic -I- potential) energy E and the mean-square value 
of its surface velocity V^, i.e., E = I^V^, and consequently 



\^{m,ijOr)f dm . 



(4) 



Here mb denotes the mass interior to the bottom boundary 
of the envelope, and i?^ and represent respectively the 
radius and the mass of the star. In practice, we normalize 
the eigenfunction ^(m, Wr) such that 



1. 



(5) 



In Fig. 6 we compare numerically computed damping 
rates with estimates from expression (2) with the right- 
hand side multiplied (arbitrarily) by 1/6. The factor 1/6 
can be obtained approximately by taking into account the 
adiabatic exponents (71 and 73) in the derivation of ex- 
pression (2), assuming a fully ionized gas. For the mode 
inertia I^^ (cf. Eq. 4) the calculations assumed adiabati- 
cally (dashed curves) and nonadiabatically (solid curves) 
computed displacement eigenfunctions ^. 

For frequencies below about half the isothermal acous- 
tical cut-off frequency uic, the results suggest a fair agree- 
ment between analytical and modelled damping rates. 
An interesting feature is the bend in the analytical solu- 
tion obtained with both the nonadiabatic and adiabatic 
eigenfunctions for the 1.45 Mq star near the frequency 
ly ~ 1.85 mHz. This property is obviously related to the 
shape of the eigenfunctions in the boundary layers of the 
convection zone, because in a polytrope the mode iner- 
tia is a smooth function of height (Gough 1995). The 
characteristic flattening of the damping rates (e.g., near 
2.8 mHz for the Sun), however, is not seen in the estimates 
from expression (2). 

6. Amplitude ratios 

A useful test of the pulsation theory, independent of 
an excitation model, is provided by comparing estimated 
intensity-velocity amplitude ratios with observations. For 
the Sun, accurate irradiance measurements exist from 
the IPHIR instrument of the PHOBOS 2 spacecraft with 
contemporaneous low-degree velocity observations in the 
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Fig. 6. Theoretical damping rates as functions of frequency 
for the Sun, a 1.25 Mq and a 1.45 ZAMS star. The curves 
are the right-hand side of expression (2) multiplied by the fac- 
tor 1/6, assuming adiabatically (dashed curves) and nonadia- 
batically (solid curves) computed mode inertia I^j. The sym- 
bols show the damping rates obtained directly from the corre- 
sponding pulsation calculations, in which were solved the fully 
nonadiabatic linearized equations using the nonlocal, time- 
dependent mixing-length model with the convection param- 
eters of Figs 4 and 5. 
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Fig. 7. Theoretical amplitude ratios for a solar model com- 
pared with observations by Schrijver et al. (1991). Computed 
results are depicted for velocity amplitudes obtained at dif- 
ferent heights above the photosphere {h = km at T = Tes) 
assuming the convection parameters used for Figs 2 and 3. The 
thick, solid curve indicates a running- mean average of the data. 

potassium line from the Birmingham instrument at Tener- 
ife (Schrijver ct al. 1991). This allows us to compare ob- 
served solar amplitude ratios with our estimated ratios 
as function of frequency. The comparison is displayed in 
Fig. 7, where the model results are depicted for veloc- 
ity amplitudes computed at different atmospheric levels. 
For moderate and high eigenfrequencies the amplitude of 
the displacement ^ increases quite steeply in the evanes- 



cent outer region of the atmosphere, where the density 
declines very rapidly. The computed velocity amplitude, 
and hence the ratio, varies by about 15% between the 
photosphere and the temperature minimum. Thus atten- 
tion has to be paid to which atmospheric level the velocity 
amplitudes are computed, i.e., at which level the displace- 
ment eigenfunctions are normalized. Observations are per- 
formed in selected Fraunhofer lines, e.g. in the neutral 
potassium line (769.9 nm) as used in the BiSON observa- 
tion (Elsworth et al. 1993), which is formed at a height 
of /i ~ 200 km above the point where the temperature 
is equal to the effective temperature (assuming the T-r 
relation derived from the model C atmosphere of Ver- 
nazza et al. 1981). The luminosity amplitudes have been 
computed at the outermost gridpoint and a correction fac- 
tor has been applied to account for the conversion to the 
measured irradiancc wavelength of A = 500 nm using the 
approximation of Kjeldsen & Bedding (1995). Observa- 
tions with a coherence greater than 0.7 are represented by 
different symbols denoting measurements of different de- 
gree I. The thick solid curve represents a running-mean av- 
erage, with a width of 300 /uHz, of the observational data. 
The theoretical ratios for h = 200 km (dashed curve) show 
reasonable agreement with the observations. 

7. Acoustical noise generation rate 

Acoustical radiation by turbulent multipole sources in 

the context of stellar aerodynamics has been considered 
by Unno & Kato (1962), Moore & Spiegel (1964), Unno 
(1964), Stein (1967), Goldreich & Keeley (1977b), Osaki 
(1990), Balmforth (1992b), Goldreich, Murray & Kumar 
(1994) and Musielak et al. (1994). 

In a pulsating atmosphere the full pulsation-convection 
equations must be derived from the fluid-dynamical equa- 
tions in which the fluid velocity includes both turbulence 
and pulsation. Balmforth (1992b) reviewed the theory of 
acoustical excitation in a pulsating atmosphere, and, fol- 
lowing Goldreich & Keeley (1977b), he derived the fol- 
lowing expression for the rate of energy injected into a 
mode with frequency cjr by quadrupole emission through 
the fluctuating Reynolds stresses: 

^1/2 



dr 



M, 



(6) 



p£QUQToS{m, LOr) dm , 



where io,Uo,To are respectively the length, velocity and 
correlation time scales of the most energetic eddies, deter- 
mined by the mixing-length model. The function iS(m, u^) 
accounts for the turbulent spectrum, which approximately 
describes contributions from eddies with different sizes to 
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the noise generation rate Pq, and which we implemented 
as did Balmforth (1992b): 

S{m, Wr) = / exp[-uj^T^/{2Ku^)^] dK , (7) 
Jo 

where k = k/ko, = u{k)/uo, k is the wavenumber of 
an eddy with velocity u{k), and fco is the wavenumber at 
the peak of the spectrum. For the computation of u{k), a 
turbulent spectrum according to Spiegel (1962) has been 
chosen. 

The emission of acoustical radiation by turbulent mul- 
tipole sources depends critically on the convective veloc- 
ity u. In homogeneous isotropic and non-decaying tur- 
bulence, acoustic emission by the fluctuating Reynolds 
stresses (quadrupolc omission) scales with the fifth power 
of the turbulent Mach number Aft = u/c (Lighthill- 
Proudman formula). Inhomogeneity and anisotropy effects 
in the overturning layers of stars give rise to monopole and 
dipole emission manifested in the fluctuation of the en- 
tropy (e.g. Goldreich & Kumar 1990). Stein & Nordlund 
(1991) and Goldreich, Murray & Kumar (1994) suggest 
that the monopole and dipole source may be as impor- 
tant as quadrupole radiation. However, previous work has 
demonstrated that the prescription is capable of roughly 
reproducing solar measurements, and so, partly for want 
of a serious theory, we stick with the expressions (6) and 
(7) here. 

8. Amplitudes 

8.1. Velocity amplitudes 

With the computations of the damping rate, r], and 
noise generation rate. Pq, the root-mean-square velocity 
at a particular level in the atmosphere may be written as 

- <«' 

The form of the turbulent spectrum S{m,u}r), given 
by Eq. (7), has a substantial effect on the predicted mode 
amplitudes. In this paper we multiplied the rhs of Eq. (7) 
by the factor 6.55 for all amplitude predictions. This em- 
pirical correction, which can be attributed perhaps to un- 
certainties in our expressions for the quadrupolc emission, 
leads to theoretical solar velocity amplitudes that have the 
same maximum value of 26.6 cm s~^ as that observed by 
the BiSON group (ChapHn et al. 1998). The results of the 
scaled theoretical mean amplitude values for a solar model 
are displayed in Fig. 8 together with the BiSON data. 

8.1.1. Main-sequence stars 

The mean velocity as a function of frequency, computed at 
a height h = 200 km above the photosphere of an evolving 
IMq star, is displayed in Fig. 9. The oscillation ampli- 
tudes become larger with age for low and intermediate fre- 
quencies, exhibiting a maximum value of Vg ~45cms~^ 




1 2 3 4 5 

u [mHz] 



Fig. 8. Velocity amplitudes for the Sun as a function of fre- 
quency. The computed values (continuous curve) are depicted 
at the photospheric level h — 200 km. The turbulent spectrum 
S{m,LOr), given by Eq. (7), has been multiplied by the factor 
6.55 to fit the maximum value of the velocity data (filled cir- 
cles) from the BiSON observations (Chaplin et al. 1998). The 
data are from a 32-month almost continuous sequence collected 
between May 1994 and January 1997, i.e., at or near the solar- 
cycle 22/23 activity minimum. The computations assumed the 
convection parameters a = 1.8, = 6^ = 300. 




u [mHz] 

Fig. 9. Velocity amplitudes for an evolving 1 Mq star as func- 
tion of frequency, depicted at the photospheric level h = 
200 km. The results are displayed for the model masses and 
convection parameters of Fig. 3. The thick curve indicates the 
results for the Sun. 

at the end of the hydrogen core-burning phase. The in- 
crease comes about because the ratio Pq/Zoj increases with 
age at the frequency of maximum mode energy, whereas 
the damping rates decrease with age at this frequency (see 
Fig. 3). 

In Fig. 10 the amplitudes are depicted for an evolving 
1.45 star, also computed at the height h = 200 km. 
For models before the characteristic 'hook' (i.e., at ages < 
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Fig. 10. Velocity amplitudes for an evolving 1.45 M© star as 
function of frequency, computed at a height h = 200 km. The 
dashed curve displays the result for the 2.44 Gy model applying 
a median filter on the amplitudes with a width corresponding 
to nine radial modes. The amplitudes are portrayed for the 
model ages and convection parameters of Fig. 4. 

2.36 Gy) in the evolutionary track (see Fig. 13) the ampli- 
tudes increase only moderately with age, mainly because 
of the increasing mode inertia at the frequency of max- 
imum mode energy and the consequent decrease of the 
ratio Pq//cj. For models older than '-^2.36Gy, the lumi- 
nosity increases fairly rapidly, leading to a steep increase 
in the turbulent Mach number and hence in the noise gen- 
eration rate Pq, and consequently mode amplitudes. Near 
the end of the hydrogen core-burning phase the theory pre- 
dicts maximum values of ~ 330 cm s~^ for the velocity am- 
plitudes. In general these maximum values coincide with 
the sharp depression in the damping rates (see Fig. 4). 
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Fig. 11. Luminosity amplitudes for ZAMS model as function 
of frequency, displayed at the outermost meshpoint of the mod- 
els. The computations assumed the convection parameters and 
model masses of Fig. 5. 



8.2. Luminosity amplitudes 

For radial modes, the imaginary part of the nonadiabatic 
displacement eigenfunctions is very small relative to the 
real part. The differences in the velocity amplitudes when 
using the adiabatic instead of the nonadiabatic displace- 
ment eigenfunctions are negligible relative to the uncer- 
tainties inherent in modelling the theory of stochastic ex- 
citation. For the estimation of the luminosity amplitudes, 
however, nonadiabatic eigenfunctions of the relative lumi- 
nosity fluctuations, SL/Lq, have to be taken into account. 
The relative luminosity amplitudes are related linearly to 
the velocity amplitudes. 
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Fig. 12. Maximum values of turbulent Mach number Mt, 

turbulent pressure fraction pt/p and convectivc growth rate 
a = 2w/i as function of model mass along the ZAMS. The com- 
putations assumed the convection parameters of Figs 4 and 5. 



8.2.1. Zero-age main-sequence stars 

The luminosity amplitudes of moderate-mass stars along 
the ZAMS are depicted in Fig. 11. The amplitudes in- 
crease monotonically with M for stars with M < 1.45 Mq, 
up to a maximum value of ~ 50ppm. For models with 
M > 1.6 Mq, amplitudes of stochastically excited modes 
decrease with M; for a 2 Mq ZAMS star the maximum 
amplitude is ~ 0.06 ppm. The dependence of the ampli- 
tude variations upon mass, or upon luminosity, may be ex- 
plained principally by the strong dependence of the acous- 
tic noise generation rate on the turbulent Mach number 
Mt. The dependence of the maximum values of the turbu- 
lent Mach number Mt and the ratio of turbulent pressure 
to total pressure pt/p upon model mass along the ZAMS 
is illustrated in Fig. 12. The computations predict the 
largest turbulent Mach numbers for models with a mass 
of ~ 1.6 Mq. The 2 Mq ZAMS star exhibits two very thin 
convection zones in the outer part of the envelope, and 
the theory predicts a maximum turbulent Mach number 
Mt < 0.1. Furthermore, the opacity and consequently the 
convective heat flux decrease with M for M > 1.6 Mq. 
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Fig. 13. Unstable modes and mean velocity amplitudes of stochastically excited oscillations. Amplitudes, evaluated at a height 
/i=200km arc depicted as contours (solid curves) labelled at the amplitude values 19, 27, 35, 45, 60, 80, 110, 150, 250, 350 
cms~^. The dotted curves are evolutionary tracks. The Sun, indicated by its symbol 0, exhibits a mean (rms) velocity of 
20.0 cms~^. Calculations have been carried out till the end of hydrogen core-burning, giving the low-temperature extremities 
of the contours. The location of the instability strips for the n = 1 and n = 2 radial modes are indicated by solid and dashed 
straight lines, respectively. 



Also indicated in Fig. 12 is the maximum value of the 
convective growth rate a = 2'w/i, scaled in units of cyclic 
frequency v. The ratio a/v influences the shape of the 
cigcnfunctions in such a way as to cause a local depres- 
sion in the damping rates rj considered as functions of v 
(cf. Gough 1997). The maximum value of a is roughly 
equal to the frequency of the local minimum of r]. 

8.3. Dependence on stellar parameters 

For those models on the main sequence for which all the 
modes are predicted to be stable, the computed maxi- 
mum velocity amplitudes of stochastically excited modes 
(evaluated at the height h = 200 km) are shown as con- 
tours on the HR diagram in Fig. 13. The 191 models 
(indicated by the diamond symbols) were generated by 
specifying the mass, luminosity and efl[ective temperature 
provided from full evolution sequences, as obtained by 



Christensen-Dalsgaard (1993). The same convective pa- 
rameters as those in Figs 4 and 5 were adopted. For more 
massive stars the maximum amplitudes exhibit peaks in 
their frequency spectrum due to the sharp dip in their 
damping rates (see Figs 4 and 10). We moderated these 
peaks by applying a median fllter to the amplitudes of all 
models with a width in frequency corresponding to nine 
radial modes (as illustrated in Fig. 10 by the dashed curve 
for a 1.45 Mq star with an age of 2.44 Gy). 

The low-temperature extremities of the contours indi- 
cate where the model hydrogen mass fractions reach 
10~^ in the core; no calculations were carried out at lower 
temperatures. The amplitudes increase steeply with lu- 
minosity, particularly for stars with mass M > 1.4 Mq, 
owing to the increase in the convective velocities with M. 
The largest amplitudes are predicted for a 1.6 Mq model 
of spectral type F2, which has a maximum velocity am- 
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Fig. 14. Luminosity amplitudes (top) and amplitude ratios 
(bottom) as function of effective temperature and model lumi- 
nosity. The amplitude ratios are displayed for velocities at two 
different atmospheric levels: the thin curves denote the results 
at h = km and the thick curves at the height h = 200 km. The 
luminosity amplitudes are computed at the outermost mesh- 
point of the models. The computations assumed model param- 
eters as given in Fig. 4. The value for the Sun (3.4 ppm) is 
indicated by its symbol. 

plitude of ~ 15 times larger than that found for the Sun. 
For this model the turbulent Mach number is also pre- 
dicted to be largest (see Fig. 12). For more massive stars 
the computations predict overstable modes (see section 9). 
Amplitudes of such overstable modes are limited by non- 
linear processes and can therefore not be estimated with 
the linear computations adopted in this paper; their val- 
ues could be much larger than the amplitude values of the 
stable stochastically driven modes considercxl here. 

For stars with L < 2L©, the velocity amplitudes of 
stochastically excited modes depend strongly on the model 
luminosity and are only weakly dependent on the effective 
temperature. The same trend is also seen in Fig. 14 for the 
luminosity amplitudes (top panel). In the lower panel of 
Fig. 14 the ratios between luminosity and velocity ampli- 
tudes, AL/AVg, are displayed, where the luminosity fluc- 
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Fig. 15. Velocity amplitudes as function of light-to-mass ra- 
tio for stochastically excited oscillations in 191 models (filled 
circles) calculated at a height /i=200km above the photo- 
sphere with the same convection parameters as for Fig. 4. The 
dashed curve indicates Kjeldsen & Bedding's scaling law (9) 
with s = 1. The amplitudes are displayed relative to the value 
found in the Sun. 

tuations, AL, arc computed at the outermost mcshpoint 
of the models. At the photospheric level {h = km) the 
amplitude ratios appear to be quite insensitive to lumi- 
nosity, and depend mainly on effective temperature. 

Based on the model results of Christensen-Dalsgaard 
& Frandsen (1983b), Kjeldsen & Bedding (1995) proposed 
a scaling relationship for solar-type velocity amplitudes as 
a function of parameters used in stellar-evolution theory. 
In particular, they proposed the scaling law 

^^(M^X (9) 

Vq \M/Mq) ' ^ ' 

with s = 1, suggesting that the velocity amplitudes scale 
directly with the light-to-mass ratio L/M oi the star. In 
Fig. 15 the velocity amplitudes versus the light-to-niass 
ratio are displayed at a height h = 200 km above the pho- 
tosphere for model calculations assuming the convection 
parameters of Fig. 4. There is a fair agreement between 
the computed amplitudes (filled circles) and Kjeldsen & 
Bedding's proposed relation (dashed line) for L/M < 3. 
For higher values of L/M the estimated amplitudes are 
predicted to be larger than Kjeldsen & Bedding's linear re- 
lation, particularly for models with masses M/M© > 1.4. 
Moreover, for these models the computed amplitudes be- 
come progressively more dependent on the model's ef- 
fective temperature and less dependent on L/M as they 
evolve along their evolutionary tracks (see also Fig. 13). 
Applying a linear polynomial fit to the estimated ampli- 
tudes in Fig. 15 suggests for the exponent s in the scaling 
law (9) a value of 1.29. At the photospheric height the 
computations suggest a value of s = 1.47. 

We should point out, however, that the convective ve- 
locities found in our models are large. As already indicated 
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Fig. 16. Luminosity amplitudes (computed at outermost mesh- 
point) versus effective temperature for models with constant lu- 
minosity. The results are displayed for different mixing-length 
parameters a (top) and metallicities (bottom). The compu- 
tations assumed the nonlocal convection parameters of Fig 4. 
The line styles are as defined in Fig. 14. Top: the thick curves 
display the results for models computed with q = 1.8 and the 
thin curves depict the amplitudes obtained with a = 2.0. In 
both model sequences the value for the metallicity Z was cho- 
sen to be 0.02. Bottom: the thick curves depict the amplitudes 
from model calculations using Z = 0.04 and the thin curves for 
Z = 0.02, assuming a = 2.0. 

in Fig. 12, the turbulent Mach number Mt becomes rel- 
atively large for models with M/Mq > 1.4. Relative to 
a local convection model, the nonlocal formulation used 
here reduces the convective velocities, although they still 
remain large. This reduction results in part from the av- 
eraging of the supcradiabatic temperature gradient over 
the eddies, which spreads the influence of this gradient's 
sharp peak in the hydrogen ionization zone over a larger 
region. 
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Fig. 17. Velocity amplitudes for an evolving l.OA/,:;, and l.SM© 
star versus the model's effective temperature. The computa- 
tions assumed the nonlocal convection parameters of Fig. 4 
and results are displayed at a photospheric level h = 200 km. 
Top: amplitudes are depicted for three values of of the mixing- 
length parameter a assuming Z = 0.02 in the computations. 
Bottom: results are plotted for three values of metallicity Z 
assuming a = 2.0 in the model calculations. 

8.4- Dependence on mixing length and metallicity 

The dependence of the luminosity amplitudes on mixing- 
length parameter a and metallicity Z is illustrated in 
Fig. 16 over a range of effective temperature for models 
with constant luminosity. The dependence of the veloc- 
ity amplitudes are illustrated in Fig. 17 for two evolving 
models with mass I.OMq and 1.3Mq. 

Increasing the mixing-length parameter results in an 
increase of both the luminosity and velocity amplitudes 
(top panels). This comes about because increasing a re- 
sults in higher convective velocities u and thus in a larger 
acoustic generation rate Pq (sec Eq. 8). Moreover, it ap- 
pears that the amplitudes become less dependent upon a 
with increasing luminosity, which might be explained by 
the decrease in convective efHcacy with model mass. 
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A similar behaviour of the amphtudes is obtained when 
the metalUcity is increased (bottom panels in Figs 16 and 
17). A larger value for the heavy-element abundance Z 
results in a higher opacity k and consequently in a larger 
convective heat flux in the upper boundary layer of the 
convection zone. Therefore the turbulent Mach number 
Mt becomes larger, and thus also do the amplitudes. 

9. Overstable modes 

For stars with logTeff > 3.85 lying more or less in the 
b Scuti instability strip, the model calculations predict 
modes to be overstable, irrespectively of adjustments to 
the convection parameters a and h. The 5 Scuti stars are 
variables with spectral types A and F in the lower part 
of the classical Cepheid instability strip, which are in the 
very interesting evolutionary phase of the main sequence 
near the end of central hydrogen burning. It was first 
shown by Zhevakin (1953) and Cox & Whitney (1958) 
that the excitation mechanism in Cepheids, which are 
core-helium-burning radial pulsators having large ampli- 
tudes of the fundamental mode (and in a few cases also 
the first overtone), is plausibly due to the opacity mech- 
anism acting in the Hell ionization zone (see also Baker 
& Kippenhahn 1962). The same mechanism is believed 
to be responsible for the excitation in 5 Scuti stars (e.g. 
Dziembowski 1995 and references therein). The oscillation 
specta of many 5 Scuti stars, however, are far more com- 
plex, involving both radial and nonradial modes with low 
amplitudes, lying often in a narrow frequency range. This 
complicates mode identification substantially (e.g. Man- 
gcney et al. 1991). The cooler ^ Scuti stars have substan- 
tial outer convection zones. Thus in these layers the pulsa- 
tionally induced fluctuations of the turbulent fluxes may 
become important for the selection mechanism of modes 
with observable amplitudes. 

The theoretically predicted order n of unstable p 
modes in sequences of evolving models of 5 Scuti stars 
are depicted by different symbols in Fig. 13 (e.g. circles 
indicate the location of models in the HR diagram for 
which the radial fundamental mode was found to be over- 
stable). The models have only a few excited modes lying 
in a narrow frequency interval, and some of them display 
radial orders in a nonconsecutive sequence (however, see 
also Houdek & Gough 1998; Michel et al. 1999). More- 
over, with increasing effective temperature the overstable 
modes shift to higher frequencies. The blue edge of the in- 
stability domain is found to shift to higher effective tem- 
peratures with increasing order; this result is consistent 
with previous model calculations (e.g. Stellingwerf 1979, 
1980; Dziembowski 1995) mainly because for these models 
convection is unimportant. Through the inclusion of the 
turbulent flux perturbations in the stability analyses the 
computations predict well defined red edges, a result which 
was previously reported by Baker & Gough (1979) for RR 



Lyrae stars. In particular, the fluctuating Reynolds stress 
(5pt is found to be the decisive contributor to the damping 
rates and thus for the return to stability at the red edge 
for low-order modes in b Scuti stars (Houdek 1997) . Only 
with the inclusion of Jpt in the c;oniputations are all modes 
found to be stable for models with effective temperatures 
satisfying log Teg < 3.85. 

10. Conclusion 

It is evident that one of the greatest deficiencies in mod- 
elling oscillations in stars with surface convection zones is 
the lack of a proper theory of convection in a pulsating 
environment. Although several attempts have been made 
in recent years to address this problem (for a review, see 
Baker 1987) none of the proposed prescriptions are any- 
thing more than phenomenological. Impressive progress 
has been made on hydrodynamical simulations of convec- 
tion, including also the interaction with pulsations (e.g. 
Stein & Nordlund 1991; Bogdan, Cattaneo & Malagoli 
1993; Nordlund & Stein 1998; Stein & Nordlund 1998). In 
particular, the work by Stein and Nordlund, including a 
realistic treatment of the physics of the outer parts of the 
convection zone, has confirmed the earlier conclusion that 
the solar oscillations are likely to be intrinsically stable 
(e.g. Gough 1980; Kumar & Goldreich 1989; Balmforth 
1992a); also, the simulations yielded estimates, similar to 
the observationally determined values, of the energy in- 
put to the modes from the stochastic driving by convec- 
tion. However, such simulations are evidently extremely 
time consuming and have so far been made in sufficient 
detail only for solar parameters. Here we have estimated 
the pulsational properties of main-sequence stars over a 
broad range of parameters, by means of a time-dependent 
non-local version of mixing-length theory. 

Perhaps the most important conclusion drawn from 
this survey is that, as in the case of the Sun, oscillations in 
solar-like stars are intrinsically damped and stochastically 
driven by convection. We note that the issue may still not 
be entirely settled, however. In particular, Cheng & Xiong 
(1997) reported calculations using Xiong's (1989) nonlocal 
formulation of mixing-length theory which predict over- 
stable solar oscillations. In fact, their computations sug- 
gest that the momentum flux perturbations destabilize all 
p modes, in complete disagreement to the results reported 
here. This discrepancy evidently deserves investigation. 
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